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Abstract 


The prompt emission mechanism of gamma-ray bursts (GRBs) is still unclear, and the time-resolved spectral 
analysis of GRBs is a powerful tool for studying their underlying physical processes. We performed a detailed 
time-resolved spectral analysis of 78 bright long GRB samples detected by Fermi/Gamma-ray Burst Monitor. A 
total of 1490 spectra were obtained and their properties were studied using a typical Band-shape model. First, the 
parameter distributions of the time-resolved spectrum are given as follows: the low-energy spectral index 
а ~ — 0.72, high- energy rds index 8 ~ — 2.42, the peak energy Е, ~ 221.69 keV, and the energy flux 
F ~ 7.49 x 10 8 erg стг . More than 8046 of the bursts exhibit the Вані low-energy spectral index oq 
exceeding the synchrotron limit (-2/3). Second, the evolution patterns of o and E, were statistically analyzed. The 
results show that for multi-pulse GRBs the intensity-tracking pattern is more common than the hard-to-soft pattern 
in the evolution of both E, and a. The hard-to-soft pattern is generally shown in single-pulse GRBs or in the initial 
pulse of multi-pulse GRBs. Finally, we found a significant positive correlation between F and Е,, with half of the 
samples exhibiting a positive correlation between F and a. We discussed the spectral evolution of different 
radiation models. The diversity of spectral evolution patterns indicates that there may be more than one radiation 
mechanism occurring in the GRB radiation process, including photospheric radiation and synchrotron radiation. 
However, it may also involve only one radiation mechanism, but more complicated physical details need to be 
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1. Introduction 


Gamma-ray bursts (GRBs) are the most intense explosions in 
the universe, which have been extensively studied since their 
discovery in the 1960s. However, there are still several 
unidentified issues during the prompt emission phase (Zhang 
et al. 2011; Kumar & Zhang 2015). For example, the 
composition of the jet (baryonic matter or Poynting flux), the 
energy dissipation and particle acceleration mechanisms 
(internal shocks or magnetic reconnection, see Zhang & 
Yan 2011), the radiation mechanisms (synchrotron radiation 
(Mészáros et al. 1994; Ravasio et al. 2018), inverse Compton 
scattering. (Kumar & McMahon 2008), and photosphere 
radiation (Meng et al. 2022). The lightcurves of GRB exhibit 
irregularity and diversity, while the spectra tend to have a few 
simpler shapes compared to the lightcurves, which can be well 
described by several experimental functions. The most typical 
of them are the Band function (Band, see Band et al. 1993), 
cutoff power law (CPL), single power law function (PL), and 
smoothly broken power law function. The above spectral 
models are usually considered non-thermal spectra. In some 
cases, the spectra of GRBs are better fitted by superimposing 
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multiple spectral functions rather than by a single spectra 
function, such as an additional power-law component (Abdo 
et al. 2009; Ackermann et al. 2010) or a thermal component 
(Ghirlanda et al. 2003; Ryde 2004, 2005; Ryde & Pe'er 2009; 
Guiriec et al. 2011; Axelsson et al. 2012; Li 2019; Zhao et al. 
2022) on the basis of a Band spectrum. Ryde (2005) fitted a 
Planck function superposition power law for some GRBs and 
found that the thermal component dominates. Abdo et al. 
(2009) fitted the spectrum of GRB 090902B with a Band 
function plus power law. Ryde et al. (2010) claimed that such a 
burst can also be fitted with a multicolor blackbody function 
model superimposed on a power law function. The thermal 
component is located on the left shoulder of the peak energy 
(below Е,) and the power-law component can be extended to 
the high and sub-low energy bands. 

The possible radiation mechanisms of the prompt emission 
phase of GRBs may consist of synchrotron radiation of 
relativistic electrons, quasi-thermal radiation from the photo- 
sphere, as well as the Comptonization of thermal photons 
(Pe'er et al. 2006; Zhang et al. 2011). Observationally, most 
of the GRB spectra behave as non-thermal spectra, so the 
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synchrotron mechanism is the preferred model to explain the 
majority of the prompt emission (Lloyd & Petrosian 2000; 
Baring & Braby 2004; Burgess 2014; Zhang et al. 2016). The 
spectral shape can be described by the Band function and the 
CPL function. In general, the low-energy spectral slope is an 
indicator to determine whether the observed radiation can be 
correlated with any of the proposed theoretical spectral 
models. The typical low-energy spectral index o allowed by 
the internal shock synchrotron model ranges from —2/3 
(slow-cooling) to —3/2 (fast-cooling). Indeed, such a value is 
much softer than many observations, living the fast-cooling 
problem of synchrotron radiation (Preece et al. 2000). 
Although the synchrotron emission in decaying magnetic 
fields could alleviate the fast-cooling problem (Uhm & 
Zhang 2014), the predicted o is still smaller than —2/3. 
The limit value —2/3 is usually thought as the “death line” of 
synchrotron radiation. Spectra with a larger than —2/3 are 
often considered as emission from the photospheres. How- 
ever, this criterion can hardly be strict. More than one 
complex emission model produce a wider range of o values. 
Models such as the subphotospheric dissipation (Rees & 
Mészáros 2005) and geometrical broadening in photospheric 
emission (Pe'er 2008) could efficiently produce a broader 
spectra. Not only that, the GRBs with —2/3 <a < 1 can be 
adequately explained with the time-dependent cooling of 
synchrotron electrons model as well (Burgess et al. 2020). 
Moreover the “fast cooling problem" can also be reconciled in 
the Internal-collision-induced Magnetic Reconnection and 
Turbulence (ICMART) model due to the turbulence heating 
(Zhang & Yan 2011; Shao & Gao 2022). Therefore, it is 
difficult to identify the emission model only by a. 

The Band function is characterized by three parameters: the 
low-energy index a, the high-energy index 6, and the peak 
energy (Ep). In previous works based on large samples, the 
typical values of the time-integrated spectrum are œ ~ — 1.0, 
B ^» — 2.0, and Ep ~ 300 keV (Preece et al. 2000; Nava et al. 
2011; Zhang et al. 2011; Goldstein et al. 2012). While for the 
time-resolved spectrum, the low-energy index is much harder 
(ағ-- 0.8, see Kaneko et al. 2006; Yu et al. 2016, 2019; Li 
et al. 2021). For bright bursts, the time-resolved spectral 
analysis provides more clues about the GRB prompt emission. 
It has been found that the characteristic parameters of the 
spectrum are not fixed, but evolve over time (Golenetskii et al. 
1983; Crider et al. 1997; Kaneko et al. 2006; Peng et al. 2009), 
making the spectral evolution pattern become another impor- 
tant indicator for studying the radiation mechanism of GRBs. 
The evolution characteristics of the time-resolved spectra E, 
and а have been extensively studied in the early days. During 
the Compton Gamma-Ray Observatory observation period, E, 
showed several different evolution patterns: (1) “hard-to-soft” 
(Norris et al. 1986); (1) “intensity-tracking” (Golenetskii et al. 
1983); (ш) “soft-to-hard” or chaotic evolution. With the launch 
of the Fermi Gamma-ray Space Telescope (Fermi) in 2008, the 
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quality of the spectral data has been improved, identifying that 
the first two patterns dominant (Lu et al. 2012). The physical 
origin of these evolved patterns has not been resolved yet. For 
single-pulse GRBs, the “hard-to-soft” pattern accounts for 
about two-thirds and the “intensity-tracking” pattern for about 
one-third (Lu et al. 2012; Yu et al. 2016, 2019). The evolution 
pattern of multi-pulse GRBs is more complex, showing hard- 
to-soft followed by intensity-tracking (first pulse “hard-to-soft” 
followed by “тасКіпе”) and tracking (including the first pulse) 
is more common (Lu et al. 2012). Recently, Li et al. (2021) 
found that E, showed a tracking pattern in 60% of their 
multipulse GRB samples. Analysis of the Fermi-LLE GRB 
spectrum revealed that the E, “tracking” behavior accounted 
for 75% of the analyzed samples (Duan & Wang 2020). As for 
the о evolution, Crider et al. (1997) first pointed out that it 
evolves over time, instead of staying the same. The o evolution 
is more complex and shows multiple evolution patterns. 
Therefore, there are relatively few studies and physical 
explanations of the a evolution. In one study, Е, and o were 
found to show a "double tracking" pattern (Li et al. 2019). 

In order to further investigate the time-resolved spectral 
properties of bright long GRBs in the prompt emission phase, 
the time-resolved spectra of 78 bright long GRBs detected by 
Fermi/Gamma-ray Burst Monitor (GBM) are analyzed in this 
work. The Band function was used to fit each burst. In addition, 
their parametric distributions are given by detailed time- 
resolved spectral analysis. The evolution patterns of a and Е, 
are also statistically analyzed. The correlations of parameters, 
such as а-Ер, Е-Е, and F-a, are obtained. 


2. Sample Selection and Spectral Analysis Methods 
2.1. Sample Selection 


Our samples come from the Fermi Gamma-ray Space 
Telescope, which was launched in 2008 June and carries two 
instruments. One is the GBM (see Meegan et al. 2009), which 
observes the energy range of 8 keV—40 MeV, and its main task 
is to search for GRBs in the universe. The other is the Large 
Area Telescope (see Atwood et al. 2009), which has an energy 
range of 30 МеУ-300 GeV. The GBM Public GRB burst 
catalog" provides up-to-date information on all triggers 
classified as GRBs since the beginning of the mission. 
Considering the requirements of the signal-to-noise ratio 
(S/N) for time-resolved analysis, we just deal with the long 
GRBs that are bright enough to have an integral flux greater 
than 2 x 10 еге cm > over the burst duration. There are 386 
GRBs satisfying such a brightness condition before 2023 
January. Furthermore, only those GRBs, whose lightcurve 
pulses are significantly identified and countable, are selected. 
Finally, 78 bright long GRBs are selected as our sample. 


^ https: / /heasarc.gsfc.nasa.gov /W3Browse /fermi /fermigbrst.html 
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2.2. Detector, Source and Background Selection 


There are 14 detectors on board Fermi/GBM, of which 12 
Sodium Iodide (Nal, n0-nb) are used to detect photons in the 
energy range of 8keV-1 MeV, and 2 Bismuth Germanate 
(BGO, bO and b1) are used to detect the photons ranged in 
200 keV—40 MeV, respectively. The Nal detectors are dis- 
tributed around the periphery of the spacecraft, with three 
detectors on each face to achieve full sky coverage. The BGO 
detectors are located on both sides of the spacecraft, which can 
similarly cover the entire sky. In general, Nal detectors п0-п5 
correspond to BGO detector bO, and Nal detectors n6-nb 
correspond to BGO detector bl. We selected data from three 
Nal detectors and one BGO detector at the optimal observation 
location. 

The data of GBM are stored in three types: the CTIME data, 
the CSPEC data, and the Time-Tagged Event (TTE) data. The 
CTIME data has a time resolution of 0.064 s with eight energy 
channels, and the CSPEC data has a time resolution of 1.024 s 
with 128 energy channels. The TTE data has the same 128 
energy channels as the CSPEC data and has a time resolution of 
2 us, containing both time and energy information for 
individual photons. Both CSPEC and TTE data have higher 
spectral resolution and can be used for spectral analysis. 
Because the TTE data has the highest temporal resolution, it is 
more suitable for time-resolved spectral analysis. The TTE and 
response files of the selected detector are used for spectral 
analysis in this work. 

The GRB source intervals we choose are generally slightly 
longer than Тоо reported in the NASA/HEASARC database, 
as this would include all relevant features of the light curve. 
Meanwhile, we select the photons in the interval of tens of 
seconds before the trigger time and the interval of tens of 
seconds after the end of the emission to estimate the 
background. A polynomial of order 0-4 is applied to fit the 
background photon counts for all 128 energy channels and the 
optimal polynomial has been determined by a likelihood ratio 
test. The polynomial is then interpolated into the source 
interval to obtain the background photon count estimate. We 
have jointly fitted the spectral data of three Nal detectors and 
one BGO detector. The spectral energy ranges of the Nal 
detectors were set to 10-30 and 40-900 keV (excluding the 
k-edge? at 33.17 keV), and the spectral energy ranges of the 
BGO detector was set to 250 keV-40 MeV. 


2.3. Lightcurve Time Binning 


The key premise for time-resolved spectrum analysis is to 
appropriately divide the lightcurve into time bins. An 
inappropriate division method may deviate from the real 
situation and thus lead to fuzzy results. There are three 
commonly used ways to divide time bins based on different 


5 https: / /fermi.gsfc.nasa.gov /ssc /data /analysis /caveats.html 
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considerations: uniform time bins, constant S/N per bin, and 
Bayesian Blocks (BBlocks, see Scargle et al. 2013). To 
minimize the variation in radiation over a time interval, 
which would mask the true spectral shape, BBlocks are 
thought to be the most appropriate method (Burgess 2014). 
We first use the BBlocks method to divide the TTE lightcurve 
of the brightest NaI detector into time bins, with setting a 
false alarm probability of ро = 0.01 (Yu et al. 2019; Li et al. 
2021). Then, the time bins are transferred and used to other 
detectors. 

The time bins divided by BBlocks have different S/Ns. 
Sometimes the S/N in a certain time interval might be too low 
to satisfy the requirements of the statistical significance fitting. 
The statistical significance S currently used is a test statistic 
that includes information on the S/N, which is applicable 
to Poisson sources with a Gaussian background (Vianello 
et al. 2018). It is necessary to choose the time bins with S > 20 
to ensure that there are enough photons to make the 
model parameters converge well (Ryde et al. 2019; Yu et al. 
2019; Li et al. 2021). In addition, there are at least four-time 
bins in each burst satisfying S > 20, so it is meaningful to 
study the spectral evolution with sufficient time-resolved 
spectra. 

The sample yielded a total of 1814 time-resolved spectra, of 
which 1490 spectra met the statistical significance of S 7 20. 
The basic information of our sample is listed in Table 1, 
including the Fermi/GBM burst ID (column 1), the detector 
used (column 2), the duration 759 (column 3), the selected 
source and background intervals (columns 4-6), the number of 
BBlocks time bins (column 7), and the number of pulses per 
burst (column 8). 


2.4. Spectral Fitting 


The typical model used in this work to fit the GRB spectra is 
the Band function (Band et al. 1993) which is described as 
follows: 


Ngana (E) -А 


а 
E | E/E, 
Еру 


(а- 8) 8 
[s = — GE , (a — B)E < E, 
E 


(a — p)Eo 2 E 


E piv piv 


(1) 


where А is the normalization factor in unit of 
photon cm ? keV! s~!; Еру is the pivot energy fixed at 
100 keV; o and 8 are the low-energy and high-energy photon 
spectral indices, respectively; E; is break energy of the 
spectrum. The peak energy E, of the vF, spectrum is related to 


Eo through Ep = (2 + a)Eo. 
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Table 1 


Basic Information of GRB Sample 
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GRB 


081009140 
081125496 
081215784 
081221681 
081224887 
090719063 
090820027 
090902462 
090926181 
091127976 
100324172 
100707032 
100719989 
101123952 
101126198 
110301214 
110625881 
110721200 
110920546 
111220486 
120119170 
120328268 
120711115 
120728434 
120919309 
130518580 
130606497 
130704560 
131014215 
140206275 
140213807 
140329295 
141028455 
150127589 
150201574 
150330828 
150403913 
150902733 
160113398 
160530667 
160905471 
160910722 
170405777 
170522657 
170808936 
170826819 
171210493 
171227000 
180113418 
180305393 
180720598 
180728728 
181227262 
190114873 
190530430 
190720613 
190727846 
190731943 


Detectors 


(n3),n4,n7,b1 
n9, (na),nb,b1 
n9, (na),nb,b1 
nO,nl, (n2),b0 
n6,n7, (n9),b1 
n6,n7, (n8),b1 
nl, (n2),n5,b0 
nO, (n1),n2,b0 
n3,n6, (n7),b1 
(n6),n9,na,b1 
nl, (n2),n5,b0 
n4,n7, (n8),b1 
n3, (n4),n5,b0 
n9, (na),nb,b1 
(n7),n8,nb,b1 
n7, (n8),nb,b1 
n7,n8, (nb),b1 
(n6),n7,n9,b1 
(n0),n1,n3,bO 
n0, (n1),n2,b0 
n7,n9, (nb),b1 
n7,n9, (nb),b1 
(n2),n8,na,b1 
nl, (n2),n5,b0 
(n1),n2,n5,b0 
(n3),n6,n7,b1 
n7, (n8),nb,b1 
n3, (n4),n5,b0 
n9,na, (nb),b1 
nO, (n1),n3,b0 
n0, (n1),n2,b0 
n7, (n8),nb,b1 
(n6),n7,n9,b1 
n7,n8, (nb),bl 
(n3),n4,n7,b0 
nl, (n2),n5,b0 
nO, (n3),n4,bO 
nO,nl, (n3),b0 
n7,n8, (nb),bl 
nl, (n2),n5,b0 
(n6),n7,n9,b1 
п1,12, (n5),b0 
n6, (n7),n9,b1 
(n0),n1,n2,bO 
п1,13, (n5),b0 
n9,na, (nb),b1 
nO, (n1),n2,bO 
n2,n4, (n5),b0 
nl, (n2),n9,b0 
nl, (n2),na,b0 
n6, (n7),n8,b1 
n3,n6, (n7),b1 
nl, (n2),na,b1 
n3, (n4),n8,b0 
nO, (n1),n2,b0 
nO, (n1),n2,b0 
nO,nl, (n3),b0 
n6,n7, (n9),b1 


Too 
(s) 
41.345 
9.280 
5.568 
29.697 
16.448 
11.392 
12.416 
19.328 
13.760 
8.701 
17.920 
81.793 
21.824 
103.938 
43.837 
5.693 
26.881 
21.822 
160.771 
39.041 
55.297 
29.697 
44.033 
100.481 
21.248 
48.577 
52.225 
6.400 
3.200 
146.690 
18.624 
21.248 
31.489 
60.929 
15.616 
153.859 
22.272 
13.568 
24.576 
9.024 
33.537 
24.320 
78.593 
7.424 
17.664 
11.008 
143.107 
37.633 
24.576 
13.056 
48.897 
6.400 
13.184 
116.354 
18.432 
6.144 
35.073 
15.872 


AT 
(s) 
0.0 to 55.0 
0.0 to 12.0 
0.0 to 10.0 
0.0 to 40.0 
0.0 to 25.0 
—1.0 to 25.0 
29.0 to 50.0 
0.0 to 30.0 
0.0 to 20.0 
—1.0 to 10.0 
0.0 to 22.0 
0.0 to 83.0 
—1.0 to 25.0 
35.0 to 160.0 
0.0 to 55.0 
—1.0 to 10.0 
0.0 to 35.0 
0.0 to 25.0 
—1.0 to 170.0 
—8.0 to 40.0 
0.0 to 70.0 
0.0 to 50.0 
0.0 to 120.0 
0.0 to 120.0 
—2.0 to 35.0 
5.0 to 65.0 
0.0 to 63.0 
—2.0 to 13.0 
1.0 to 6.0 
0.0 to 50.0 
0.0 to 20.0 
—1.0 to 30.0 
0.0 to 50.0 
0.0 to 75.0 
0.0 to 25.0 
0.0 to 170.0 
0.0 to 45.0 
0.0 to 25.0 
20.0 to 60.0 
—2.0 to 20.0 
0.0 to 50.0 
0.0 to 40.0 
0.0 to 88.0 
0.0 to 10.0 
0.0 to 30.0 
0.0 to 15.0 
0.0 to 175.0 
0.0 to 55.0 
2.0 to 32.0 
0.0 to 20.0 
0.0 to 60.0 
9.0 to 21.0 
0.0 to 20.0 
0.0 to 120.0 
0.0 to 25.0 
—1.0 to 9.0 
0.0 to 41.0 
0.0 to 22.0 


AT x91 
(s) 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to—10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 
—25.0 to —10.0 


АТые2 
(s) 
60.0 to 80.0 
20.0 to 40.0 
15.0 to 30.0 
50.0 to 65.0 
30.0 to 60.0 
35.0 to 50.0 
65.0 to 80.0 
35.0 to 50.0 
30.0 to 45.0 
20.0 to 40.0 
30.0 to 45.0 
100.0 to 115.0 
40.0 to 60.0 
175.0 to 200.0 
60.0 to 75.0 
20.0 to 40.0 
50.0 to 65.0 
35.0 to 50.0 
175.0 to 190.0 
65.0 to 80.0 
80.0 to 95.0 
775.0 to 90.0 
145.0 to 160.0 
200.0 to 215.0 
45.0 to 60.0 
85.0 to 100.0 
75.0 to 95.0 
25.0 to 40.0 
20.0 to 40.0 
70.0 to 90.0 
35.0 to 50.0 
40.0 to 65.0 
60.0 to 85.0 
100.0 to 125.0 
50.0 to 75.0 
225.0 to 240.0 
50.0 to 75.0 
50.0 to 75.0 
75.0 to 100.0 
40.0 to 60.0 
75.0 to 100.0 
60.0 to 85.0 
100.0 to 115.0 
15.0 to 30.0 
50.0 to 65.0 
25.0 to 40.0 
185.0 to 200.0 
60.0 to 75.0 
40.0 to 55.0 
35.0 to 50.0 
100.0 to 115.0 
35.0 to 50.0 
35.0 to 50.0 
130.0 to 145.0 
35.0 to 50.0 
15.0 to 35.0 
50.0 to 75.0 
35.0 to 50.0 


Spectra(Ns>20) 
(М) 
28(20) 
9(6) 
26(22) 
16(14) 
10(7) 
14(11) 
21(18) 
51(47) 
26(24) 
22(17) 
11(7) 
16(12) 
19(12) 
48(33) 
10(9) 
18(14) 
384) 
9(8) 
10) 
35(20) 
14(10) 
23(20) 
33(25) 
50(41) 
9(5) 
20(16) 
43(36) 
19(16) 
29(27) 
29(17) 
13(10) 
26(18) 
14(11) 
23(17) 
25(25) 
41(34) 
16(12) 
16(12) 
15(10) 
22(19) 
17(11) 
20(18) 
21(14) 
10(6) 
39(37) 
15(11) 
14(9) 
43(40) 
22021) 
12(10) 
75(70) 
22(18) 
25(22) 
49(45) 
54(54) 
13(9) 
30(19) 
18(14) 
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Table 1 
(Continued) 
GRB Detectors Too АТ. AT ke,1 AT Spectra(Ns>20) Pulse 
(8) (5) (s) (s) (М) (М) 
200101861 п1,13, (п5),50 9.984 0.0 to 20.0 —25.0 to —10.0 30.0 to 45.0 25(19) 2 
200125864 nO,nl, (n2),bO 5.824 —1.0 to 10.0 —25.0 to —10.0 20.0 to 40.0 33(30) 5 
200313071 n4,n7, (n8),b1 13.568 —1.0 to 15.0 —25.0 to —10.0 25.0 to 40.0 16(12) 1 
200412381 n6, (n7),n8,b1 6.080 0.0 to 14.0 —25.0 to —10.0 25.0 to 40.0 30(25) 2 
200826923 n1,n2, (n5),bO 7.424 0.0 to 11.0 —25.0 to —10.0 20.0 to 35.0 18(15) 1 
200829582. 14,16, (n8),b1 6.912 15.0 to 30.0 —25.0 to —10.0 40.0 to 55.0 18(15) 1 
201016019 n3,n4, (п5) ЬО 2.944 0.0 to 11.0 —25.0 to —10.0 20.0 to 35.0 27(18) 1 
201216963 n9, (na),nb,b1 29.953 0.0 to 40.0 —25.0 to —10.0 55.0 to 70.0 33(31) 2 
210406949 п1,13, (n5),bO 19.712 0.0 to 30.0 —25.0 to —10.0 45.0 to 60.0 17(16) 3 
210610827 (n9),na,nb,b1 55.041 10.0 to 90.0 —25.0 to —10.0 120.0 to 135.0 17(14) 3 
210714331 n9,na, (nb),bl 42.369 0.0 to 45.0 —25.0 to —10.0 55.0 to 70.0 15(11) 1 
210801581 n9, (na),nb,b1 13.824 —1.0 to 17.0 —25.0 to —10.0 25.0 to 40.0 14(10) 2 
211019250 nO, (n1),n2,bO 47.361 0.0 to 52.0 —25.0 to —10.0 60.0 to 75.0 15(10) 1 
220304228 n8, (na),nb,b1 31.489 0.0 to 40.0 —25.0 to —10.0 55.0 to 70.0 11(8) 1 
220426285 nO,nl, (n2),bO 5.632 0.0 to 10.0 —25.0 to —10.0 20.0 to 35.0 24(23) 2 
220527387 n6, (n7),n8,b1 10.496 0.0 to 15.0 —25.0 to —10.0 25.0 to 40.0 23(20) 4 
220910242 17,18, (nb),b1 4.224 0.0 to 11.0 —25.0 to —10.0 20.0 to 35.0 23(18) 3 
221022955 n3, (n4),n5,b0 31.744 0.0 to 50.0 —25.0 to —10.0 60.0 to 75.0 20(17) 2 
221023862 nO, (n1),n2,bO 39.169 0.0 to 50.0 —25.0 to —10.0 60.0 to 75.0 32(26) 1 
221209243 nl, (n2),n5,b0 4.160 0.0 to 8.0 —25.0 to —10.0 15.0 to 35.0 15(11) 1 


Note. Fermi/GBM burst ID (column 1), utilized detector with the brightest response indicated in parentheses (column 2), duration represented by Тоо (column 3), 
source intervals AT; (column 4), and background intervals АТ, з, before the trigger time (column 5), and background intervals АТь, 2 after the end of emission 
(column 6), number of BBlock time bins (column 7), number of time bins with statistical significance greater than or equal to 20 enclosed in parentheses within 


column 7, and count of pulses observed in each burst (column 8). 


The Multi-Mission Maximum Likelihood Framework pack- 
age (3ML,° Vianello et al. 2015) is used for spectral fitting and 
parameter estimation in this work. There are two approaches 
available in 3ML for fitting data and models: the maximum 
likelihood estimation (MLE) method and the Bayesian method. 
In this study, we employ Bayesian parameter estimation, which 
involves accessing the posterior distribution of model para- 
meters using a specific sampling algorithm. The optimal model 
parameters are estimated from the posterior probability 
distribution, which is evaluated by combining the prior 
distribution with the likelihood function that quantifies how 
well the model matches the observed data. Typical spectral 
parameters from the previous Fermi/GBM catalog are selected 
as priors, where the normalization factor (A) is a logarithmic 
uniform distribution (log U), the low and high energy indices а 
and 8 are Gaussian (С), and the peak energy E, of the vF, 
spectrum is a logarithmic normal distribution (log N). We apply 
emcee (Foreman-Mackey et al. 2013) in this work for posterior 
sampling, which is a Markov Chain Monte Carlo (MCMC) 
sampling algorithm. We specify the number of walkers to 20, 
and the number of global samples to 10,000. During MCMC 
sampling, it requires a certain number of samples for the 
Markov chain to reach convergence and achieve a steady-state 


6 https:/ /threeml.readthedocs.io/en /stable /index.html 


of the parameter distribution. Therefore, the samples that did 
not reach the steady-state distribution in the previous period 
should be discarded. In this study, we excluded the initial 2596 
of 10,000 MCMC samples for each parameter sampling and 
only considered the last 75% as representative samples from 
the posterior distribution that reached convergence. 

There are two main aspects of spectral fitting. The first one is 
parameter estimation, which aims to determine the best-fitting 
parameters and their uncertainties for a given model. The 
second aspect involves evaluating the concordance between the 
model and the data, commonly referred to as the goodness-of- 
fit (GOF). Traditionally, the GOF is evaluated using a reduced 
chi-square (x7), denoted as x^ /dof, where dof represents the 
degree of freedom. A good fit is usually indicated by a 
value close to 1. However, this approach may lack reliability 
and accuracy when the data deviates from a normal distribu- 
tion. In the case of MLE or Bayesian methods, alternative 
statistical measures can be used to assess the GOF of the 
data. For example, Cstat (Cash 1979) and Pgstat, where 
Pgstat = —21nL(0), with L(0) representing the likelihood 
value as a function of the free parameter 0. The minimum of the 
—|InL(0) function corresponds to the maximum likelihood. 
Different statistics are utilized in Fermi data analysis due to the 


7 https:/ /heasarc.gsfc.nasa.gov /xanadu /xspec /manual /XSappendixStatistics.html 
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small counts in each time bin, which renders it inappropriate to 
treat each bin as a single observation from a normal 
distribution. Instead, a Poisson distribution is more suitable. 
It is widely acknowledged that as the bin counts increase 
sufficiently, the Poisson distribution approximates a normal 
distribution. At this point, Pgstat is expected to closely 
approximate X^, and it can be anticipated that —2 In L(0)/dof 
will approach 1. However, if the bin counts are insufficiently 
large, there is no basis for expecting —2 In L(0)/dof equal to 1. 


3. Spectral Analysis Results 


We employed the Bayesian approach (emcee) to conduct 
spectral fiting on 1490 time-resolved spectra with a 
significance of S > 20 obtained from 78 bright long GRBs. 
As an illustrative example, Table 2 presents the spectral fitting 
results for GRB 171227000. Similar to the method employed 
by Li (2019), we utilized Pgstat as the statistical measure 
instead of x7. Consequently, —21nL(0)/dof does not 
typically converge to unity. As previously elucidated, this 
occurrence is not uncommon due to insufficient photon counts 
in each bin. 


3.1. Parameter Distribution 


The global parameter distribution, including a, 8, E, and the 
derived parameter flux (F), F is the best-fit model integrand in 
the energy range of 1 keV—40 MeV, is depicted in Figure 1. The 
Gaussian profiles were employed to fit the a and ( distributions 
(N= u + о, where u represents the average value and с 
denotes the corresponding standard deviation). On the other 
hand, due to their approximate lognormal, a lognormal profile 
was used to fit the distributions of E, and Ғ. The overall 
sample exhibited mean values and standard deviations as follows: 
a = —0.72 + 032, 8—-—242 + 0.39, logj(Ep/keV) = 
log,o(221.69) + 0.41, log, )(F /(erg cm? s^!) = log) 9(7.49e — 
6) + 0.59. These findings are consistent with the results of 
previous studies (Poolakkil et al. 2021). 

We conducted an analysis of 1490 spectra, where 669 
spectra (44.9%) of them displayed o values surpassing the 
synchrotron radiation limit (a > — 2/3). According to Acuner 
et al. (2020), spectra with a > — 0.5 indicate a preference for 
the photosphere model. Therefore, based on their criteria, we 
identified 410 spectra (27.596) that exhibit a stronger inclina- 
tion toward photosphere emission. 

Subsequently, the hardest low-energy index @max was 
identified for each burst, following a distribution as depicted in 
Figure 2, oq, = —0.34 + 0.35. Notably, 67 bursts (86%) 
exhibited an Qa, value surpassing the synchrotron radiation 
limit, which significantly exceeded the proportion observed in 
the overall time-resolved spectrum. According to the criteria of 
Acuner et al. (2020), 56 bursts (71.7%) have Qmax > —0.5,a 
pronounced inclination toward photosphere emission within 
this spectral range. 
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The presence of a harder low-energy index suggests that the 
simple synchrotron model is inadequate to fit the spectrum, and 
it might be usually considered as an indicator of the 
photosphere model. If a harder o is owed to the photosphere 
emission, a thermal component can be expected in the 
spectrum. However, as referred in Section 1, an inference of 
this kind could not so indubitably be drawn, considering the 
different details of the radiation process (e.g., Lundman et al. 
2013; Burgess et al. 2020). Indeed, it might be misleading to 
judge the radiation mechanism with relying solely on whether 
the low-energy index exceeds the synchrotron radiation limit. 
A more reasonable method is to fit the GRB data directly using 
the physical models rather than empirical functions and 
compare the goodness of fit (e.g., Meng et al. 2018). For 
example, by directly fitting the radiation model involving time- 
dependent cooling of synchrotron electrons to the observed 
data of GRBs, Burgess et al. (2020) pointed out that the 
synchrotron emission call still be suitable for many GRBs, 
despite their a larger than —2/3. Some other works also prove 
the necessity and feasibility of implementing the directly fitting 
to data with physical model (Zhang et al. 2016; Yang et al. 
2023), but it is beyond the scope of this work. 


3.2. Spectral Evolution 


Previous studies have primarily focused on the evolution of 
parameters, particularly the peak energy Ер, which has been 
extensively investigated in early research. The evolution of E, 
exhibits various patterns: (1) the hard-to-soft pattern, (ii) the 
intensity-tracking pattern, and (iii) the soft-to-hard pattern or 
chaotic evolution. Subsequent studies have highlighted the 
prevalence of the first two patterns. In the case of multi-pulse 
GRBs, however, spectral evolution patterns become more 
intricate, and recent research has revealed a higher occurrence 
of rough tracking rather than smooth tracking (Duan & 
Wang 2020; Li et al. 2021). Moreover, there is also noticeable 
evolution behavior in the low-energy index o, which is more 
complex compared to the evolution of Ep. Li et al. (2021) found 
that half of their samples exhibited tracking-patterns of а 
evolution. On the other hand, the high-energy index 3 does not 
display a clear population trend, and its greater chaos-like 
behavior has not been statistically analyzed before. 

The evolution of spectral parameters, such as a, p, Е, and 
VF, flux for 6 GRBs as representative examples, is illustrated in 
Figure 3. Based on the observed trends along with the 
lightcurve, we categorized the evolution patterns of o and Ep 
for all samples. Specifically, the patterns were classified into 
different categories: hard-to-soft (h.t.s), intensity-tracking (i.t), 
hard-to-soft followed by intensity-tracking (h.t.s to 1.0, and 
other patterns. In this classification, “hard” and "soft" denote 
larger and smaller values of о and Е parameters, respectively. 
The detailed results can be found in columns 3 and 4 of 
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Table 2 
Time-resolved Spectral Fit Results of GRB 171227000 

tobh S а 8 Ер Е BIC/-In(posterior) dof 
(s) ~ (s) (keV) (erg cm ? 57!) 

0.0 ~ 6.348 22.23 —0.16+01 2.361013 467.321516 3.072056 x 1076 5197.8 /2586.63 463 
6.348 ~ 14.39 46.21 —0.43*008 —2.81:022 411.0321 7 3.37033 x 10 5 5590.1/2782.75 463 
14.39 ~ 15.099 23.39 0461009 2.421015 562.4417633 8.861146 x 10-9 2232.5/1103.95 438 
15.099 ~ 15.548 282 —0.49 008 242575046 627.34212:9 13.6511 4! x 107% 1783.2/879.3 398 
15.548 ~ 16.821 68.64 —0.59 1008 ҮЗІП 1097.97+8972 32.89+145 x 1076 3400.1/1687.77 463 
16.821 ~ 17.341 59.64 0.621006 215042 1284.91 +436 51.514303 x 107° 2197.9 /1086.64 413 
17.341 ~ 17.648 36.2 —0.66 +0.05 —2.481011 891.03 1133! 28.3713 o x 1076 1461.1 /718.25 375 
17.648 ~ 17.82 44.69 -0.46104 —2.34+0.08 995.8107? 71.337322 x 107% 1020.7/498.05 313 
17.82 ~ 18.185 47.03 —0.47 005 -2.37109 660.9555; 34.5413 i x 10-5 1736.0/855.72 383 
18.185 ~ 18.647 83.44 —0.670 03 -2.391007 1446.55110774 114.671416 x 107% 2298.4/1136.91 400 
18.647 ~ 20.433 170.86 0.61109! 5955004 1249.664508 97.96*13, x 107% 4366.1/2170.77 463 
20.433 ~ 20.836 74.95 -0.64:09 -2,52109 868.011-23 58.3573:48 x 1076 2006.2/990.81 387 
20.836 -- 21.554 73.28 —0.66 18:03 —2.68 +911 799.447393] 33.314138 x 1076 2643.8 /1309.6 438 
21.554 ~ 22.212 81.8 —0.65 003 —2.731059 854.0 203 41.6424) x 107% 2589.3/1282.35 436 
22.212 ~ 22.728 65.1 -0.56 103 —2.26100% 515.6 48544 36.1728 x 1076 2157.6/1066.51 413 
22.728 ~ 22.985 35.72 —0.65*007 -2/71047 502.48:542! 15.3511 x 1076 1152.8/564.13 354 
22.985 -- 24.215 55.76 —0.65 +0.04 —2.35*04 433.86 135, 12.534137 x 1075 3101.5/1538.44 463 
24.215 ~ 25.597 4221 —0.627008 -2.61:047 372432188 6.22107 x 107% 3217.4/1596.39 463 
25.597 ~ 25.942 31.04 —0.69+006 —2.52+015 562.03:26-66 13.29*79, x 1076 1525.4/750.43 380 
25.942 ~ 26.489 27.53 —0.75+8!! —2.027$4, 376.37+1998 9.993242 x 107% 1881.6/928.5 419 
26.489 ~ 26.738 40.58 0.711007 2.61023 827.99116228 29.847333 x 107% 1270.9/623.15 351 
26.738 ~ 27.252 40.04 —0.72 +0.06 -2.69102, 521.315222 12.71414 x 10 9 1957.4/966.39 413 
27.252 ~ 27.403 31.79 0.541007 22109 448125192 24.59148, x 107% 743.7/359.56 292 
27.403 -- 27.557 52.74 —0.63 +0: o —2287*015 901.06 *12415 81.07+648 x 1076 998.2/486.78 294 
27.557 ~ 28.019 59.26 —0.71+005 —2.68*043 554.21 * $206 23.381224 x 10 6 1967.8/971.6 402 
28.019 ~ 28.712 63.39 —0.77+004 —2.43+012 506.6372;% 19.374791] x 1076 2469.4/1222.41 437 
28.712 ~ 29.102 37.0 0.771008 2.56102, 452.4113837 12.27413 х 1076 1621.9/798.65 385 
29.102 ~ 29.373 45.47 0.692002 —2.62+018 516.55+5531 22.5114 x 1076 1380.0/677.73 360 
29.373 ~ 29.91 43.72 —0.727045 —2.36+01, 340.4322:06 11.37*?3 x 1076 1981.9/978.66 418 
29.91 ~ 31.369 40.43 -0.591049 29143008 201.05101, 5.831092 x 1076 3138.9/1557.15 463 
31.369 ~ 31.757 28.09 —0.79+011 2.46102, 307.1515088 6.19138 х 10 6 1501.2/738.28 384 
31.757 ~ 32.395 221 —0.67*013 -22670H 19611461 3.381106 x 107% 1976.9/976.17 436 
32.395 ~ 38.267 42.74 —0.867005 -2.43 5043 134.6 1593 1.51023 x 1076 5008.1/2491.77 463 
38.511 ~ 39.062 36.6 —0.83^005 -2.431043 201.8713998 5.371112 x 1076 1887.5 /931.47 421 
39.062 ~ 40.274 44.06 КҮЗЕТІ -2.07109 149.42120. 5.861193 x 107% 2890.1/1432.77 463 
40.274 ~ 42.928 39.29 —0.93 19:14 -2,1710%9 121.26 181 2.445044 x 1076 3886.1/1930.74 463 
42.928 ~ 45.357 26.46 -0.861022 -2.25109 91.847145] 1.41026 x 10-6 3710.2/1842.81 463 
45.357 ~ 46.238 21.02 —0.94* 011 —2415* 031 117.4413 1.727972 x 107° 2390.7/1183.07 449 
46.238 ~ 49.915 31.2 —1,04828 —2.63*021 83.84* 3950, 0.94*018 x 1076 4302.0/2138.72 463 
49.915 ~ 51.705 29.18 —1.0+019 07712019 78.41726 1.081043 x 1076 3268.7/1622.03 463 


Note. Тһе start and stop times of the BBlock time bins (column 1), the significance S of time bins (column 2), the optimal parameters for the Band model (columns 
3-5), the derived energy flux (column 6), the Bayesian information criterion BIC and the -In(posterior) values (consistent with the —In L (0) values) (column 7), as well 


as degrees of freedom (column 8). 


Table 3. The overall statistics for the evolution patterns are 
summarized in Figure 4 and Table 4. 

Regarding the evolution of Ер among our total sample of 78 
long GRBs, the intensity-tracking pattern was exhibited by 35 
bursts (45%), the hard-to-soft followed by intensity-tracking 
pattern was displayed by 23 bursts (2996), the hard-to-soft 


pattern was shown by 17 bursts (22%), while three bursts (4%) 
exhibited other evolution patterns. The dominant evolution 
patterns for E, were the hard-to-soft pattern and the intensity- 
tracking pattern. Interestingly, in some GRBs both of these 
patterns coexisted, resulting in the hard-to-soft followed by 
intensity-tracking pattern. 
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Figure 1. The distributions of the time-resolved spectral parameters a (upper left panel), 8 (upper right panel), E, (lower left panel), and F (lower right panel). The 
low-energy index (o) and high-energy index () are fitted with the Gaussian profile, and the peak energy (Ер) and flux (F) are fitted using a log-normal distribution. 


The red line represents the fitted line of the parameter distribution histogram. 
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Figure 2. The distribution of the hardest low-energy index (Qmax) in each GRB, 
and the red line represents a Gaussian fit to the Qmax distribution. 


In terms of the evolution of o, the observed patterns were 
characterized by a high degree of complexity and unpredict- 
ability. The majority of our sample showed tracking patterns, 
hard-to-soft patterns, and hard-to-soft followed by intensity- 


tracking patterns. The evolution of o was roughly assessed, 
with 36 bursts (46%) displaying a rough intensity-tracking 
pattern, 14 bursts (1896) showing the hard-to-soft evolution 
pattern, 17 bursts (2296) exhibiting the hard-to-soft followed by 
intensity-tracking pattern, and 11 GRBs (14%) in the sample 
demonstrating other evolution patterns. Specifically, we 
observed several cases with different evolution patterns, such 
as the soft-to-hard followed by intensity-tracking (e.g., GRB 
1402134), the intensity-tracking followed by soft-to-hard (e.g., 
GRB 200125B), and one flat evolution pattern (GRB 110920). 

It is worth noting that the underlying physics processes 
governing single-pulse and multi-pulse bursts may differ 
significantly, thus it is crucial to investigate the spectral evolution 
patterns separately for these bursts. In our sample data set, there 
аге 29 single-pulse GRBs. For E, evolutionary behavior within 
this subset, three main patterns emerged 16 bursts (5596) 
exhibited the hard-to-soft evolution, eight bursts (2896) showed 
the intensity-tracking evolution, and four bursts (14%) displayed 
the hard-to-soft followed by intensity-tracking evolution pattern. 
One burst (GRB 110721A) demonstrated a  hard-to-soft 
evolution followed by soft-to-hard pattern. In regards to the 
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Table 3 
Evolution of Parameters E, and o and Parameter Relations 

GRB a — 2/3 E, a а-Е, Е-Е, Е-о Pulse 

(Evolution) (Evolution) Type(r) Type(r) Type(r) (N) 
081009140 not all it it 2p(0.5) 2p(0.94) 1(0.45) 2 
081125496 yes hts it 1(0.38) 2p(0.94) 2p(0.46) 1 
081215784 not all h.ts to i.t it 1(0.37) 2p(0.9) 2p(0.64) 3 
081221681 not all h.t.s to i.t h.t.s to i.t 2p(0.79) 2p(0.92) 1(0.71) 2 
081224887 not all hits hts 2p(0.94) 2p(1) 2p(0.94) 1 
090719063 not all h.ts to i.t hts 2p(0.85) 2p(0.95) 2p(0.78) 1 
090820027 not all it hts 2p(0.63) 2p(0.9) 1(0.73) 1 
090902462 not all it h.ts to it 3(0.01) 2p(0.87) 3(—0.01) 3 
090926181 not all it it 3(0.32) 2p(0.76) 3(0.28) 2 
091127976 not all h.ts to i.t it 1(0.53) 2p(0.39) 2p(0.45) 3 
100324172 not all h.t.s to i.t h.t.s 2p(0.64) 2p(0.82) 1(0.57) 2 
100707032 not all hits it 2p(0.95) 2p(0.97) 2p(0.99) 1 
100719989 not all it h.ts to it 1(0.28) 2p(0.8) 2p(0.48) 3 
101123952 not all h.ts to it to h.ts it 3(0.48) 2p(0.82) 2p(0.64) 5 
101126198 no it it to s.t.h 2p(0.55) 2р(0.89) 1(0.54) 1 
110301214 not all h.t.s to i.t h.t.s to i.t 2p(0.64) 2p(0.75) 1(0.55) 2 
110625881 not all h.t.s to i.t it 2p(0.53) 2p(0.79) 3(0.46) 3 
110721200 no h.t.s to s.t.h it 1(0.02) 1(0.52) 1(0.67) 1 
110920546 yes hts flat 2n(—0.7) 2p(0.98) 3(—0.62) 1 
111220486 no h.t.s to i.t it 3(0.19) 2p(0.65) 3(0.37) 2 
120119170 no h.ts to i.t h.ts to it 2p(0.68) 2p(0.76) 2p(0.5) 1 
120328268 not all h.ts to i.t it 2p(0.64) 2p(0.85) 2p(0.75) 2 
120711115 no it h.ts to it 3(—0.15) 3(0.21) 2p(0.67) 2 
120728434 not all it c 2n(—0.92) 2p(0.79) 2n(—0.61) 5 
120919309 not all it h.ts to it 2n(0.67) 2p(0.9) 2p(0.5) 1 
130518580 not all it h.t.s to i.t 1(0.37) 2p(0.91) 3(0.46) 1 
130606497 not all it h.ts to it 3(-0.41) 2p(0.87) 3(—0.14) 4 
130704560 not all hts it 2n(0.89) 1(0.73) 2p(0.85) 3 
131014215 not all it hts 2n(0.44) 2p(0.75) 1(0.23) 2 
140206275 not all h.ts to it h.ts to it 2n(0.83) 1(0.88) 2p(0.89) 2 
140213807 not all it s.t.h to it 3(0.09) 2p(0.85) 1(0.21) 2 
140329295 not all h.t.s to i.t it 3(0.2) 2p(0.78) 1(0.44) 2 
141028455 not all hts i.t to s.th 1(0.33) 1(0.58) 2p(0.64) 1 
150127589 not all it to hts h.t.s 2p(0.91) 2p(0.57) 1(0.57) 2 
150201574 not all it it 2p(0.7) 2p(0.89) 2p(0.9) 2 
150330828 not all it h.ts to it 2p(0.66) 2p(0.65) 1(0.49) 5 
150403913 no h.ts to i.t it 2p(0.46) 1(0.60) 2p(0.93) 2 
150902733 not all h.ts to i.t it 1(0.30) 2p(0.63) 2p(0.66) 1 
160113398 yes hts it 1(0.51) 1(0.81) 2p(0.84) 1 
160530667 not all it it 2p(0.7) 2p(0.89) 2p(0.8) 1 
160905471 not all it it 2p(0.64) 2p(0.92) 2р(0.71) 2 
160910722 not all hits it 2p(0.79) 1(0.49) 2p(0.70) 1 
170405777 not all it i.t to hts 3(0.29) 2p(0.66) 3(0.48) 3 
170522657 yes h.ts to it hts 1(0.43) 1(0.26) 2n(—0.54) 2 
170808936 not all it it 1(0.81) 2p(0.88) 2p(0.82) 3 
170826819 not all it h.t.s 2p(0.43) 2p(0.84) 2n(0.55) 3 
171210493 not all h.t.s s.t.h to i.t 1(-0.25) 2ра) 3(-0.25) 1 
171227000 not all it h.ts to it 2p(0.63) 2(0.93) 2p(0.53) 3 
180113418 not all h.ts to i.t flat to s.th 2n(—0.19) 1(0.69) 2n(—0.34) 2 
180305393 yes h.t.s h.t.s to i.t 2n(0.71) 2p(0.77) 3(0.38) 1 
180720598 no h.ts to i.t it 1(0.5) 2p(0.8) 2p(0.72) 5 
180728728 no it h.t.s 3(0.16) 2p(0.91) 3(0.24) 1 
181227262 not all h.t.s to i.t h.t.s to i.t 2p(0.48) 2p(0.83) 1(0.30) 2 
190114873 not all it flat to i.t 2p(0.87) 2p0.95) 1(0.89) 3 
190530430 not all it it 3(0.4) 2p(0.81) 2p(0.77) 3 
190720613 no h.ts to i.t it 1(0.65) 2p(0.83) 2p(0.74) 3 
190727846 no it it 1(0.42) 2p(0.70) 3(0.5) 4 
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Table 3 
(Continued) 

GRB a — 2/3 Ер а а-Е, Е-Е, Е-а Pulse 

(Evolution) (Evolution) Type(r) Type(r) Type(r) (N) 
190731943 not all hits to it h.t.s to it 1(0.42) 2p(0.71) 2p(0.88) 1 
200101861 not all it it 1(0.4) 2p(0.68) 2p(0.79) 2 
200125864 not all it it to s.t.h 2p(0.5) 2p(0.88) 2p(0.61) 5 
200313071 not all h.t.s it to s.t.h 1(0.32) 2p(0.70) 1(0.75) 1 
200412381 not all it it 1(0.32) 2p(0.81) 2p(0.70) 2 
200826923 not all it it 2p(0.51) 1(0.61) 3(0.58) 1 
200829582. not all h.t.s it 1(0.55) 1(0.62) 2p(0.86) 1 
201016019 not all h.t.s it 1(0.68) 1(0.78) 2p(0.93) 1 
201216963 no it it 3(0.31) 2p(0.73) 3(0.31) 2 
210406949 no h.ts to it it 2p(0.89) 1(0.61) 2p(0.58) 3 
210610827 yes it h.t.s 2p(0.54) 2p(0.94) 3(0.42) 3 
210714331 not all h.t.s it 3(0.35) 2p(0.66) 2p(0.79) 1 
210801581 no it it 2p(0.60) 2p(0.71) 2p(0.52) 2 
211019250 not all h.t.s h.t.s 2p0.95) 2p(0.98) 2p(0.97) 1 
220304228 yes hts h.t.s 2p(0.9) 2p(0.76) 2р(0.64) 1 
220426285 yes h.t.s to i.t h.t.s to i.t 2p(0.86) 1(0.75) 1(0.63) 2 
220527387 not all it it 3(0.33) 2p(0.59) 2p(0.82) 4 
220910242 not all h.ts to it it 2p(0.49) 2p(0.87) 2p(0.60) 3 
221022955 not all it h.t.s 2p(0.44) 2p(0.74) 3(0.16) 2 
221023862 no it h.t.s to i.t 3(0.23) 2p(0.92) 2p(0.420 1 
221209243 not all h.t.s h.t.s 2p(0.89) 2p(0.85) 2p(0.89) 1 


Note. Spectral proportions of œ across synchrotron death lines (column 2), Е, evolution (column 3), a evolution (column 4), types of parameter relations and 


Spearman rank coefficients (columns 5-7) were analyzed. 


evolution of o in single-pulse GRBs, 11 bursts (3896) showed an 
intensity-tracking pattern, six bursts (21%) displayed the hard-to- 
soft followed by intensity-tracking pattern, and seven bursts 
(2496) demonstrated the hard-to-soft trend. 

Among the 49 multi-pulse GRBs, 27 bursts (5596) exhibited 
ап intensity-tracking pattern for the evolution of Ер, while 17 
bursts (35%) demonstrated the hard-to-soft followed by 
intensity-tracking pattern. One burst (GRB 130704) exhibited 
the hard-to-soft pattern, and four bursts (8%) displayed other 
evolution patterns. For the evolution of а, 25 bursts (51%) 
showed the intensity-tracking pattern. 11 bursts (2396) showed 
the hard-to-soft followed by intensity-tracking pattern, with the 
hard-to-soft pattern usually appearing in the first pulse. Seven 
bursts (1496) showed hard-to-soft pattern. These evolution 
patterns observed for both o and Е, align with previous studies 
(Yu et al. 2019; Duan & Wang 2020; Li et al. 2021). Lu et al. 
(2012) suggested that all subsequent pulses exhibiting the 
intensity-tracking behavior were attributed to superimposed 
hard-to-soft pulses; however, it is noteworthy that the intensity- 
tracking pattern was also identified within the first pulse itself, 
indicating that this phenomenon represents a genuine char- 
acteristic of E, evolution. 


3.3. Global Parameter Relations 


The correlations among parameters play a crucial role in 
understanding the characteristics of prompt radiation from 
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GRBs. Previous studies have primarily focused on investigat- 
ing relationships between various parameters. An early 
discovery was the Golenetskii correlation, which established 
a strong correlation between flux (F) and peak energy (Ep) 
(Golenetskii et al. 1983). In this study, we further examined the 
global correlations among selected sample parameters. The 
overall relationships between four groups of parameters are 
depicted in Figure 5: the F—o relation (top left panel), the F—E, 
relation (top right panel), the а-Е, relation (bottom left panel), 
and the a—( relation (bottom right panel). We quantified these 
parameter correlations using the Spearman correlation coeffi- 
cient (ғ), where r- 0.7 indicates a strong correlation, and 
r< 0.4 indicates a weak correlation. 

Our results indicate a strong positive correlation for the F-E, 
relation in the log-log plot, with a correlation coefficient of 
r—0.702(p < 107%, thus confirming the validity of the 
Golenetskii correlation. The fitted Е-Е, relationship is given as 

log F = 0.971003 


logio (Ep) — 7.37*006, (2) 


where the units of E, and Ғ are keV and erg cm ? 571, 


respectively. Additionally, we also separately analyzed the 
Е-Е, relationship for single-pulse and multi-pulse GRBs, 
shown in Figure 6, resulting in two distinct relationships: 


F = 0.971004 10g1o(Ep) — 7.457013 


g (single), 
logio (Ep) — 7.3770 


logo 
logo F = 0.997003 


(multi). (3) 
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Figure 4. Histogram of parameter evolution, left panel displays the temporal evolution of E,, while right panel illustrates the evolution of a. The first group shows 
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pattern, "i.t" signifies an intensity-tracking pattern, "h.t.s to i.t" indicates a hard-to-soft followed by intensity-tracking pattern, and "other" encompasses other 


evolution patterns. 


Table 4 

Statistical Results of Ej and o Evolution 
GRB Sample h.t.s it h.t.s to i.t Other? 

(М?) N(PCT^) N(PCT) N(PCT) N(PCT) 
E, Evolution 
Single-pulse 29 16(5596) 8(28%) 4(14%) 1(3%) 
Multi-pulse 49 1(2%) 27(55%) 17(35%) 4(8%) 
АП 78 17(22%) 35(45%) 23(29%) 3(4%) 
а Evolution 
Single-pulse 29 (2496) 11(3896) 6(21%) 5(17%) 
Multi-pulse 49 7(14%) 25(51%) 11(2396) 6(1296) 
АП 78 14(18%) 36(46%) 17(22%) 11(14%) 
Notes. 


* Hard-to-soft (h.t.s); intensity-tracking (i.t); hard-to-soft followed by intensity- 
tracking (h.t.s to i.t); other patterns (other). 

Number of samples. 
с 

Percentage. 


Their Spearman rank correlation coefficients аге 
r — 0.606(p < 10 ^) and r= 0.734(p < 10 ^, respectively. It 
can be observed that both single-pulse and multi-pulse GRBs 
exhibit power-law exponents ~1 for their respective F—E, 
relationships without any apparent difference. The absolute 
values of the correlation coefficients for the other three 
relations are below 0.4, indicating weak correlations. 


3.4. Individual Parameter Relations 


We conducted a comprehensive analysis of parameter 
relationships for each GRB, and some examples are presented 
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in Figure 3. Following the categorization scheme proposed by 
Yu et al. (2019), we classified the individual relationships into 
three categories: (i) non-monotonic relations (type 1), which 
include both positive and negative power-law segments; (11) 
monotonic relations that can be described by a single power 
law, including monotonic positive correlations (type 2p) and 
monotonic negative correlations (type 2n); and (iii) relation- 
ships exhibiting no apparent trend (type 3). We quantified 
these parameter relationships using Spearman’s rank coeffi- 
cients. The specific relation types and corresponding Spear- 
man coefficients for each burst are summarized in columns 
5-7 of Table 3. The distribution of relation types and 
correlation coefficients for individual bursts is illustrated in 
Figure 7. 


3.4.1. Е-Е, Relation 


From the results presented in Figure 3, Table 3, and Figure 7, 
it can be observed that the relationship between flux (F) and 
peak energy (Ep) in GRBs predominantly exhibits either a 
monotonic correlation (type 2) or a non-monotonic correlation 
(type 1) when plotted on a logarithmic scale. Among the 
individual Е-Е, relationships, 63 bursts (81%) are categorized 
as monotonically positively correlated (type 2p), while 14 
bursts (18%) demonstrate а non-monotonic trend (type 1), 
additionally, one burst does not exhibit apparent trend (type 3). 
Within the sample set, strong correlations (0.7 « r « 1) are 
observed in 59 cases (76%), moderate correlations 
(0.4 « r « 0.7) are found in 16 instances (20%), while three 
bursts (4%) indicate weak correlations (0 < r < 0.4). 
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Figure 6. The global F—E, relationship for single-pulse and multi-pulse GRBs, 
orange and blue data points represent single-pulse and multi-pulse GRBs, 
respectively. The fitted lines are represented by dashed and solid lines 
correspondingly. Spearman's correlation rank coefficients for single-pulse (r) 
and multi-pulse (rœ) GRBs were calculated. 


3.4.2. F-a Relation 


Regarding the F—a relation, it is observed that 42 bursts 
(5396) exhibit a monotonically increasing trend (type 2p), while 
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four bursts (596) display a monotonically decreasing trend (type 
2n). Additionally, 16 bursts (21%) demonstrate a non- 
monotonic trend (type 1), and another 16 bursts (21%) do 
not exhibit a clear trend pattern (type 3). Analyzing Spearman's 
rank coefficient reveals that out of these bursts, 29 bursts (3796) 
show a strong correlation, whereas, for 33 bursts (42%), there 
is a moderate correlation. Furthermore, in the case of nine 
bursts (12%) show correlations, seven bursts (9%) display 
negative associations. Notably, seven bursts are found to have a 
negative correlation, out of which three bursts demonstrate a 
moderately negative correlation (—0.7 « r « — 0.4) and four 
bursts show a weakly negative correlation (—0.4 « r « 0). 


3.4.3. о-Е, Relation 


Regarding the о-Е, relation, it is observed that 33 bursts 
(42%) exhibit a monotonically positive correlation when 
plotted on a linear-log scale (type 2p), while nine bursts 
(1296) display a monotonically negative correlation (type 2n), 
Additionally, 21 bursts (27%) show non-monotonic relation- 
ships (type 1), and for 15 bursts (1996), no significant trend is 
demonstrated (type 3). Analyzing the statistical results of the 
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Figure 7. Top panel: histograms illustrating the parameter relations types, 
specifically а-Е,, Е-Е, and F-a relations for groups 1, 2, and 3, respectively. 
Bottom panel: distribution histograms of Spearman's rank coefficient r for the 
parameter relations. 


correlation coefficients, it is revealed that 19 bursts (24%) 
exhibit strong correlations, 33 bursts (42%) show moderate 
correlations, 20 bursts (2696) display weak correlations, 6 
bursts (8%) demonstrate negative correlations (—1 < r < 0). 
Notably, GRB 120728B shows a strong negative correlation 
(r= — 0.92). It can be observed that the a—E, relation does not 
have a clearly dominant type of relationship, and most of the 
samples exhibit relatively weak correlations. 


4. Summary and Discussion 


In this work, we performed a detailed time-resolved spectral 
analysis of 78 bright long GRBs detected by Fermi/GBM. Our 
selected sample yielded 1490 time-resolved spectra that satisfy 
the statistical significance S > 20. The Band function was 
employed to fit all time-resolved spectra. First, we statistically 
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analyzed the parameters (a, B, Ep), the derived parameter (F) of 
all time-resolved spectra as well as the maximum (hardest) of 
the low-energy index. Second, we investigated the evolution 
patterns of both the peak energy E, and the low-energy index 
a. Finally, we examined correlations among parameter 
relations F—E,, Е-а, and а-Е,. 

The fitted results of the time-resolved spectral parameters are as 
follows: a = — 0.72 + 0.32, 8 — — 2.42 + 0.39, log,)(Ep/keV)= 
log,(221.69) + 0.41, log, (F/(erg сш 28-1) = logjg(7.49e—6) + 
0.59. The value of a іп our sample is slightly harder (—0.72) than 
the typical value of the previous works (—0.8) but remains within 
the synchrotron limit. The hardest low-energy index Qmax in each 
burst has 8696 surpasses the synchrotron limit. 

As for the distribution of spectral evolution patterns, there 
exists a slight disparity between multi-pulse GRBs and single- 
pulse GRBs. For multi-pulse GRBs, one can conclude that the 
intensity-tracking pattern is more common than the hard-to-soft 
pattern for the evolution of both the peak energy Ep and the 
low-energy index o. For single-pulse GRBs, the evolution 
pattern of peak energy Ер was more common for the hard-to- 
soft than the intensity-tracking pattern. It should be pointed out 
that the universally existed intensity-tracking pattern of E, in 
multi-pulse GRBs might be a natural aspect of the positive 
relationship of F—E, showed before. 

The parameter relationship between the global sample and 
the individual samples was also investigated, revealing a robust 
positive correlation between E, and F in the entire data set. In 
the individual samples, a significant majority (8196) of the 
F-E, relations exhibit a monotonic positive correlation, 
indicating an intrinsic association between flux and peak 
energy. Additionally, approximately half (53%) of the Е-а 
relations exhibit a monotonic positive correlation. The para- 
meter correlations have been investigated in some previous 
works, the F—E, relationship has been found to be strongly 
correlated both in the time-integrated spectra of large GRB 
samples (e.g. Golenetskii et al. 1983; Borgonovo & 
Ryde 2001) and the time-resolved ones of an individual GRB 
(e.g., Lu et al. 2012; Yu et al. 2019; Li et al. 2021). The F-« 
relationship was found to be positively correlated in some GRB 
time-resolved spectra (Ryde et al. 2019; Yu et al. 2019; Li et al. 
2021). Our results are generally consistent with the study of 
103 pulses from 38 multi-pulse GRBs by Li et al. (2021). For 
the Е-Е, relationship, they found that 74 pulses (71%) 
exhibited monotonic positive correlations, and for ће F—a 
relationship, they found that 69 pulses (6796) exhibited 
monotonic positive correlations. 

The Band component of most observed GRB spectra is 
widely believed to originate from synchrotron radiation. There 
are at least two models proposed for interpreting the prompt 
emission of GRBs: the internal shocks (IS) model, as described 
by Paczynski (1986), Goodman (1986), Shemi & Piran (1990), 
Rees & Mészáros (1992, 1994), Mészáros & Rees (1993) and 
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the ICMART model, as described by Zhang & Yan (2011). 
Both models are capable of producing synchrotron radiation. 
In the internal shock model, the peak energy Ер о 
YiL RM +z)! (Zhang & Mészáros 2002), where L is 
the “wind” luminosity of the ejecta, ye is the electron Lorentz 
factor in the emission region, R is the emission radius, and z is 
the redshift of the GRB. Despite predicting a positive 
correlation between peak energy (Ер) and luminosity (flux), 
the power law index of the theoretical model does not align 
with the observed F—E, relations. The positive correlation of 
the Е-а relation can be attributed to subphotospheric heating 
within the flow with varying entropy as explained by Ryde 
et al. (2019). According to their findings, during the peak phase 
of the peak of the light curve, a high entropy causes the 
photosphere to approach its saturation radius, leading to a 
narrow spectrum with intense emission. Conversely, during the 
pulse decay phase when entropy decreases, weaker emission 
and broader spectrum are expected as the photosphere moves 
away from its saturation radius. 

As mentioned above, the hard-to-soft pattern of Е, evolution 
is abundantly present in single-pulse GRBs and the first pulse of 
multi-pulse GRB, and the intensity-tracking pattern of Е, 
evolution is present in a fraction of single-pulse GRBs and a 
large number of multi-pulse GRBs. Both two evolution patterns 
have already been predicted by different emission models. The 
IS model can give a tracking behavior, since Ер ос p (Zhang 
& Mészáros 2002). A hard-to-soft evolution pattern of E, 
evolution is predicted by the ICMART model (Zhang & 
Yan 2011). On the other hand, the photosphere model can also 
reproduce an E,-tracking pattern (Deng & Zhang 2014). In 
partial multi-pulse GRBs, the hard-to-soft pattern and the 
intensity-tracking pattern of E, evolution coexist in different 
pulses of a same GRB, and there is a transition between these 
two evolution patterns. As for the physical interpretation, two 
different possibilities could be drawn from the coexistence and 
transition of two evolution patterns in a same GRB. One is that 
two different emission mechanisms separately produce two Е, 
evolution patterns and transfer from one to the other. Another 
point of view is that only one complicated emission model 
(which considers finer physical details) can also produce the 
coexistence and transition of two evolution patterns. Uhm et al. 
(2018) demonstrated that a synchrotron model within a bulk- 
accelerated emission region can successfully reproduce both a 
hard-to-soft evolution and an intensity-tracking pattern of Ep. 
Furthermore, Gao et al. (2021) demonstrated through numerical 
simulations that the synchrotron model can achieve a tracking 
pattern of E, in two cases, one is that the cooling process of 
electrons is dominated by adiabatic cooling or synchrotron self- 
Compton (SSC, see Derishev et al. 2001; Geng et al. 2018) 
+adiabatic cooling at the same time, the other is that the 
emitting region is accelerated and dominated by SSC cooling. 
Otherwise, а hard-to-soft pattern is normally expected. 
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Additionally, the photospheric emission model from a structured 
jet can reproduce a hard-to-soft pattern and an intensity-tracking 
pattern of E, evolution (Meng et al. 2019). In addition, Shao & 
Gao (2022) showed that the ICMART model can produce the 
hard-to-soft pattern and the intensity-tracking pattern of Ер 
evolution. For one ICMART event, its intrinsic evolution of E, 
is a hard-to-soft pattern. Considering one observed single-pulse 
could be formed by overlapping many sub-pulses produced by 
multiple ICMART events, the resultant evolution of Е can 
exhibit an intensity-tracking pattern. In multi-pulse GRB 
produced by multiple ICMART events, the two evolution 
patterns of E, can coexist with a variety of pattern transitions. 
This scenario can also reproduce most of the Е, evolution 
patterns found in this paper. Noting the majority of hard-to-soft 
pattern in sing-pulse GRBs, the ICMART seems to be a very 
competitive model. Іп summary, the diversity of spectral 
evolution patterns indicates that there may be more than one 
radiation mechanism occurring in the GRB radiation process, 
including photospheric radiation and synchrotron radiation. 
However, it may also involve only one radiation mechanism, but 
more complicated physical details need to be considered. 
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